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COSMOS  462  (1971- 106 A) : ORBIT  DETERMINATION  AND  ANALYSIS 

by 

I ^ ' Doreen  M.  C.  Walker 

K 

Summary 

Cosmos  462  (19  71  — 1 Oft A)  was  launched  on  3 December  1971  into  an  orbit  inclined 
at  65.75°  to  the  equator,  with  a perigee  height  of  230  km  and  apogee  height  of 
1800  km.  The  satellite  remained  in  orbit  for  40  months  and  decayed  on  4 April  1975. 
Orbital  parameters  have  been  determined  at  85  epochs,  using  the  RAE  orbit  refine- 
ment program,  DROP,  with  6635  radar  and  optical  observations,  including  197  from 
the  Hewitt  cameras.  The  average  standard  deviation  in  eccentricity  and  inclin-  V 

ation  corresponded  to  a positional  accuracy  of  about  100  m.  In  addition,  orbits 
of  similar  accuracy  were  determined  daily  for  the  last  15  days  of  the  life,  from 
2000  NORAD  observations. 

During  its  slow  decay,  the  orbit  passed  through  14:1,  29:2  and  15:1 
resonances  with  the  Earth's  gravitational  field.  The  variations  in  inclination 
and  eccentricity  at  these  resonances  have  been  analysed  in  detail  to  evaluate 
lumped  geopotential  harmonic  coefficients  of  order  14,  29  and  15.  v 

The  variation  of  inclination  between  resonances  has  been  analysed  to  obtain 
four  values  of  the  average  atmospheric  rotation  rate  A at  heights  of  200-250  km 
in  1972-1975.  The  values  of  A show  a seasonal  dependence,  being  greater  in 
winter  than  in  summer,  and  the  average  rotation  rate  is  lower  than  in  the  1960s, 
being  near  1.0  rev/day.  Analysis  of  the  inclination  in  the  last  15  days  of  the 
satellite's  life  indicates  a weak  west-co-east  wind  at  high  latitude  (54-62°N). 

The  variation  of  perigee  height  has  been  analysed  to  obtain  24  values  of 
density  scale  height  H , including  eight  in  the  last  15  days.  Comparison  with 
values  from  CLRA  1372,  shows  a bias  difference  of  only  1%  and  rms  difference  of 
10 Z;  so  CIRA  1972  provides  a good  approximation  to  the  values  of  H in  1972-1975. 
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Cosmos  462  entered  orbit  on  3 December  1971  and  was  one  of  a pair  of  satel- 

1-3 

lites  launched  by  the  Russians  to  test  high-speed  interception  . Cosmos  462 
was  the  hunter  satellite  and  exploded  after  passing  the  target  satellite.  Cosmos 
459.  This  explosion  occurred  3.5  hours  after  launch  and  the  largest  piece  of  the 
satellite  remaining  in  orbit  was  designated  1971-I06A. 

After  this  experiment,  1971—1 06A  had  an  orbital  period  of  about  105  minutes, 

perigee  height  230  km,  apogee  height  1800  km  and  inclination,  i , 65.7°;  the 

satellite  remained  in  orbit  for  40  months,  without  further  disturbance, and  decayed 

4 

naturally  on  1975  April  4.90  . 

Early  in  its  lifetime,  1971—1 06A  was  selected  for  high-priority  observing 
by  the  British  optical  and  radar  tracking  stations,  including  the  Hewitt  cameras 
at  Malvern  and  Edinburgh;  1971—1 06 A was  in  an  orbit  which  would  be  useful  for 
determining  the  atmospheric  rotational  speed  from  the  decrease  in  orbital  inclin- 
ation, and  air  density  at  heights  near  200  km  could  be  evaluated  from  its  decay 
rate.  i* 

The  orbit  of  1971-I06A  has  been  determined  from  all  available  observations 
with  the  aid  of  the  RAE  orbit  refinement  program  PR0PJ,  in  the  PR0P6  version^,  at 
85  epochs  during  its  40-month  life.  In  addition  15  daily  orbits  at  the  end  of 
the  life  were  determined  from  NORAD  observations.  'This  Report  describes  the  orbit 
determinations,  and  the  analyses  of  variations  in  the  orbital  parameters.  The 
changes  in  inclination  and  eccentricity  at  14:1,  29:2  and  15:1  resonances  have 
been  analysed  to  determine  lumped  geopotential  coefficients  of  order  14,  29  and 
15.  (There  was  no  detectable  perturbation  at  31:2  resonance.)  The  values  of 
inclination  between  the  resonances  have  been  analysed  to  determine  the  atmospheric 
rotation  rate  at  heights  near  200  km.  The  accurate  values  of  perigee  height 
obtained,  which  have  already  been  used  in  evaluating  the  air  density^,  have  been 
analysed  to  determine  the  atmospheric  density  scale  height.  The  residuals  of  the 
observations  have  been  used  to  assess  the  accuracy  of  each  observing  station. 

2 THE  OBSERVATIONS 

The  orbit  of  1971—1 06A  has  been  determined  at  85  epochs  from  6635  observa- 
tions. A breakdown  of  the  number  and  type  of  observations  used  on  each  of  the 
85  runs  is  given  in  Table  1 on  page  4. 

The  observations  can  be  divided  into  six  groups  (see  Table  1),  the  most 
accurate  being  those  from  the  Hewitt  cameras  at  Malvern  (M)  and  Edinburgh  (E) . 

These  observations,  which  were  available  on  43  transits,  29  from  Malvern  and  14 


089 


from  Edinburgh,  usually  have  an  accuracy  of  2 seconds  of  arc  m position  and 
1 millisecond  in  time. 


Hewitt 

camera 


5E,4M 

IM 


5E , 7M 


5M.5E 

I5E, 15M 
5M 


Table  1 

Sources  of  the  observations  used  in  each  run 


Source  of  observations 


Cape 

kinetheodol ite 


Visual 
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The  observations  in  tbi>  second  group,  made  by  the  kinetheodoli  te  at  the 
South  African  Astronomical  Observatory,  have  an  accuracy  of  I minute  of  arc. 

These  observations  are  very  valuable  because  they  can  greatly  improve  the  orbital 
accuracy,  especially  the  values  of  eccentricity,  as  they  are  made  in  the  southern 
hemisphere,  whereas  northern-hemisphere  observations  are  predominant. 

The  third  group  consists  of  visual  observations  made  by  volunteer  observers 
reporting  to  the  Radio  and  Space  Research  Station  (now  Appleton  Laboratory), 

Slough,  and  to  the  Moonwatch  Division  of  the  Smithsonian  Astrophysical  Observa- 
tory. These  visual  observations  usually  have  accuracies  between  1 and  4 minutes 
of  arc.  This  group  accounted  for  about  20Z  of  the  total  number  of  observations 
available. 

The  fourth  group  of  observations  are  those  made  by  British  radar  stations. 
About  80  observations  were  made  by  the  radar  tracker  at  RRF,  Malvern,  and  the  rest, 
around  2250,  by  the  radar  trackers  at  RAF,  Fylingdales. 

The  fifth  group  consists  of  US  Navy  observations,  supplied  by  the  US  Naval 
Research  Laboratory.  Some  2550  observations  were  available  with  a topocentric 
accuracy  of  about  2 minutes  of  arc.  The  final  group  of  120  observations  comes 
from  the  theodolite  at  Jokioinen,  Finland,  with  accuracies  of  about  5 minutes 
of  arc. 

1 THF  ORBITS  PBTA1NKP  AND  THK  PBSLKVAT 10NA1.  ACCURACY 
3. 1 The  orbits 

Orbits  were  determined  at  85  epochs  fairly  evenly  spaced  over  the  satellite's 
life,  and  the  orbital  elements  at  each  epoch  are  listed  in  Table  2 on  pages  *8-41 
with  the  standard  deviations  below  each  value.  The  epoch  for  each  orbit  is  at 
00  hours  on  the  day  indicated.  In  the  1’ROPb  model*  the  mean  anomaly  M is 
fitted  by  a polynomial  of  the  form 

2 3 4 5 

M = M.  + M t + M + Mt  + M t + M t , (1) 

0 1 2 3 4 5 

where  t is  the  time  measured  from  epoch  and  the  number  of  M-coef f icients  used 

depends  on  the  drag.  For  a high-drag  orbit  like  that  of  1971-lObA,  the  number  of 

coefficients  to  be  used  is  found  by  trial  and  error.  Best  results  were  obtained 

using  Mq  to  M,.  , the  full  complement  of  coefficients  allowed  in  the  PROF  model, 

for  32  of  the  85  orbits;  33  orbits  required  M_  to  M,  ; lb  orbits  M.  to  M, 

0 4 0 3 

and  the  remaining  4 needed  only  coefficients  to  M.,  . 
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Hu-  orbits  .ill  tit  tlio  ot"  . i vat  i .mis  m .1  satistactory  manner.  witli 
partuneter  indicating  t lie  measure  ot  lit,  ranging  between  O.V>q  ami  0.10,  where  > 
is  the  sum  ot  tlio  squares  ot  the  weighted  residuals  divided  b\  the  uumbet  ot 
degrees  ot  freedom;  and  the  weighted  icsidtial  is  detined  as  the  value  lot  tight 
ascension  ot  dee  1 i nat  unit  given  by  the  orbit,  minus  the  value  given  by  the 
observation,  divided  by  the  ; v r error  assigned  to  tin-  obsei vat  ion.  No 
part i eu tar  ditticulty  m tilting  was  expert eneed  at  the  times  ot  magnetic  slot  ms, 
when  the  drag  is  liable  to  vary  irregularis.  lot  example,  the  orbit  at  the  lime 
ot  the  major  solar  disturbance  in  early  August  I97J  lifted  well. 

Hewitt  camera  obsei  vat  ions  were  available  tot  inclusion  in  .’1  ot  the  orbits 
determined,  see  Tables  1 and  For  these  -hi  orbits  the  standard  deviation  in 

eccentricity  ranged  1 1 otu  0. 00000.’  10  0.0000  '■  with  an  average  ot  0.0000IJ,  i vni  ■ 
ponding  t o about  HO  m in  perigee  or  npo.ee  height.  The  standard  deviations  in 
inclination  v.u  ied  trout  0.0001  to  0.0011  giving  an  average  sd  ot  O.OOOh  , also 
equivalent  to  about  80  m in  distance.  The  other  el  orbits,  those  without  Hewitt 
camera  observations,  bad  standard  deviations  in  inclination  varying  from  0.000. 
to  0.00  10l  with  an  avers  0 ot  0.001  ‘ , equivalent  to  some  170  m in  distance,  more 
than  double  the  average  sd  ot  tin  orbits  with  Hewitt  camera  observations.  The 
improvement  in  accuracy  lor  eccintricitv  is  less  significant,  being  about  15*. 

The  same  inereasi  in  aeeurac>  t o r eecent 1 i.  -ity  is  achieved  on  those  orbits  using 
the  tape  k 1 lie  t heodo.l  i t < observat  i ons  . This  eentinus  the  opinion  already  exptessc 
that,  while  the  Hewitt  c.uuera  observations  greatly  improve  the  accuracy  ot  the 
orbital  inclination  land  the  otbet  orbital  parameters  to  a lesser  degreed,  the 
southern  hemisphere  observat ions  have  a great  iniluetiee  in  detining  the  shape  ot 
the  orbit,  as  shown  by  the  increased  iccurac>  in  eccentricity. 

Fig  I shows  the  values  ot  eccentricity,  from  Table  d:  the  gradual  decrease 
due  to  drag  is  the  dominant  feature;  tor  most  satellites  the  oscillation  due  to 
odd  .onal  harmonics  is  usual  tv  also  important,  but  it  is  so  small  as  t o be  impel 
ceptihle  for  1^71  lOeA,  because  the  inclination  it''.  / 1 is  close  to  the  value 

(66.l°f  where  the  effect  of  the  third  onal  harmonic  is  cancelled  bv  the  ettects 

of  the  fifth  and  higher  harmonics  (see  Ret  8).  Hie  maxima  and  minima  in  slope  ot 
the  curve  in  Fig  I coincide  with  the  minima  and  maxima  in  the  variation  ot  the 
satellite's  perigee  height,  which  occur  when  the  perigee  is  at  the  equatoi  ot  at 
maximum  latitude.  Fig  •»  ot  Kol  7 shows  that  perigee  height  has  maxima  at 

MJl)  4 I S 7 5 , 4 1 9d0  and  < 0,  and  minima  at  MJ1)  41410,  41760,  41?  105  arnd  4.’4d0. 

(The  triangles  in  Fig  1 indicate  values  trout  the  dail\  orbits  for  the  last  l1'  dav 
of  the  life,  see  section  4.) 
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Fig  2 shows  the  values  of  inclination,  i , from  Table  2 (and  Table  4). 

These  values  of  inclination  are  irregular,  and  before  any  meaningful  conclusions 
can  be  drawn  from  the  variation  of  inclination,  the  values  must  be  cleared  of  all 
perturbations  (see  sections  5 and  6).  However,  the  perturbing  effects  of  the 
14th-  and  15th-order  resonances  are  discernible  in  the  observational  values,  at 
dates  centred  on  MJD  41659  and  42302  respectively. 

The  accuracy  of  other  orbital  parameters  is  much  as  expected.  Most  of  the 
values  of  right  ascension  of  the  node,  Q , have  standard  deviations  of  0.001° 
or  0.002°,  generally  slightly  higher  than  the  sd  in  the  inclination.  In  the 
first  half  of  the  life,  the  argument  of  perigee,  w , is  accurate  to  0.01°,  but 
as  usual  the  sd  tends  to  vary  as  e ' and  increases  considerably  towards  the 
end  of  the  life  as  e ->  0 . 

The  variation  of  is  for  1971—1 06A  is  much  slower  than  for  most  satellites, 
because  the  inclination  is  quite  close  to  the  critical  value  of  63.4°,  for  which 
w = 0 . Fig  3 shows  the  values  of  w from  Table  2:  the  argument  of  perigee 
decreases  by  about  per  day,  and  perigee  does  not  quite  complete  2 revolutions 
during  the  40  months  of  the  satellite's  life. 

Nearly  all  the  values  of  M(  , are  accurate  to  better  than  I part  in  10^, 
and  consequently  nearly  all  the  values  of  semi  major  axis,  a , have  standard 
deviations  of  between  1 and  3 m. 

The  values  of  M,,  , which  provide  a direct  measure  of  the  air  drag^,  are 
mostly  accurate  to  better  than  J%:  they  are  plotted  in  Fig  4 and  have  been  fully- 
utilized  and  discussed  in  Ref  7. 

3.2  Accuracy  of  the  observations 

The  orbit  refinement  program  proceeds  by  rejecting  observations  which  do 
not  fit  well  (weighted  residuals  >3c):  a total  of  5449  observations  out  of  the 
original  6635  were  accepted  in  the  final  orbits. 

The  residuals  of  the  observations  have  been  obtained  using  the  ORES  computer 
9 

program  and  sent  to  the  observers.  The  accuracies  of  selected  observing  stations, 
with  five  or  more  observations  accepted  in  the  orbit  determination,  are  listed  in 
Table  3,  page  9,  along  with  the  number  of  accepted  observations.  The  US  Navy 
observations  from  station  29  are  geocentric,  and  if  they  were  given  in  the  same 
form  as  the  other  ( topocentr ic)  observations,  their  angular  rms  residuals  would 
increase  by  a factor  of  about  five.  The  total  rms  residuals  in  Table  3 are  much 
as  expected:  0.04  minutes  of  arc  for  the  Hewitt  cameras;  1.3  minutes  of  arc  for 
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the  Cape  kinetheodol ite;  and  from  2.0  minutes  of  arc  upwards  for  the  visual 
observers  and  theodolites. 

Table  3 

Residuals  for  observing  stations  with  more  than  five  observations 
accepted  in  the  orbit  determination 


Stat i on 


Number  ol 
obse  rva i i ons 
accepted 


Range 


Rms  residuals 

Minutes  of  arc 


1 

US  Navy 

2 

US  Navy 

3 

US  Navy 

4 

l'S  Navy 

5 

US  Navy 

6 

US  Navy 

29 

US  Navy 

81) 
2 594  ( 

Lender  zee  1 

4 1 4 

Capetown 

725 

Bucharest 

1114 

Mi  skolc 

1 9t>3 

Jokioinen 

2285 

Far nham 

2303 

Malvern 

Hewitt  camera 

2304 

Malvern  radar 

2406 

Dubl in 

2414 

Bournemou  th 

2419 

Tremadoc 

2420 

Wi  1 lowbrae 

2421 

Malvern  4 

2430 

Stevenage  4 

2437 

Warrington 

2502 

Sudbury 

2513 

Colchester 

2528 

Aldershot 

2534 

Edinburgh 
Hewitt  camera 

2539 

Dymchurch 

2550 

M as i rah 

2577 

Cape 

k inetheodol ite 

2596 

Akrotiri 

4126 

Groningen 

4130 

Denekamp 

8597 

Adel  aide 

1 .7 
1 .5 

4.7 

5.8 

10.  I 

1 .6 
0.02 


2 

.8 

2 

.8 

2 

. 1 

2 

2 

2 

.3 

O 

.5 

10 


0KB ITS  FOR  TUI'  FIFTEEN  DAYS  BEFORE  DECAY 

After  the  orbital  determinations  described  in  section  3 were  completed,  a 
further  2000  observations  were  provided  by  the  assigned  and  contributing  sensors 
of  the  North  American  Air  Defense  Command  (N0RAD)  Space  Detection  and  Tracking 
System  (ST AD ATS)  for  the  last  IS  days  before  decay.  With  the  aid  of  these 
observations  it  has  been  possible  to  determine  orbits  at  the  end  of  the  satellite's 
life  more  accurately  and  at  closer  intervals  than  lias  been  possible  in  the  past. 

Fifteen  further  orbits  were  determined  at  daily  intervals  from  1975  March 

21.0,  using  these  NORAD  observations  together  with  other  observations  previously 

used  ia  orbits  83-85  of  Table  2.  The  orbital  elements  are  listed  in  Table  4 on 

page  a 2 , with  the  standard  deviations  below  each  value.  As  before,  the  epoch 

for  each  orbit  is  at  00  hours  on  the  day  indicated.  In  II  of  these  15  orbits, 

only  the  coefficients  M , and  M,  were  required  in  the  polynomial  for  mean 

anomaly,  equation  f 1 ) . This  is  unusual,  because  the  full  set  of  coefficients  is 

generally  needed  near  decay  (as  in  Table  2):  here,  however,  the  observations  for 

each  orbit  extend  over  no  more  than  24  hours,  so  that  t ..  0.5  and  terms  in  t *, 

4 5 

t and  t have  much  less  effect  than  when  t > 1 as  for  the  orbits  of  Table  2.  i' 

For  example,  if  the  set  of  coefficients  in  orbit  84  of  Table  2 were  correct,  the 
value  of  M,t  1=  0.05t**)  would  be  <0.003C  for  t <0.5  . So  the  term  would  not 

•4 

be  significant  for  the  one-day  orbit,  and  M,  would  not  be  needed.  The  values 

*4 

of  c,  i and  M,  from  Table  ■*  have  been  added  at  the  end  of  Figs  1,  2 and  4 
respectively . 

The  last  three  orbits  in  Table  2 (orbits  83,  84  and  851  are  virtually 
independent  of  the  corresponding  orbits  in  Table  4 (orbits  E,  J and  N) , because 
the  latter  orbits  include  only  a few  of  the  original  observations,  which  are 
greatly  outweighed  by  the  large  number  of  NORAD  observations.  So  it  is  interest- 
ing to  compare  the  elements.  is  based  on  only  30  observations  spread 

over  5.1  days,  and  requires  all  six  M-coef f icients,  so  that  10  parameters  are 
being  determined  from  30  observations.  This  might  almost  be  called  a recipe  for 
unreliability.  The  corresponding  one-day  orbit,  E , with  78  observations  and 
only  M ^ , M|  and  M.,  , should  be  much  more  reliable.  Comparison  shows  that 
the  values  of  e,  0 and  M,  in  orbit  81  are  inconsistent  with  those  ol  orbit  E, 

but  the  values  of  i,  u>,  M and  Mj  are  within  the  combined  sd . So  orbit  83 

emerges  from  the  comparison  surprisingly  well.  Orbit.  S-!  is  based  on  only  lb  obser- 
vations spread  over  2.1  days,  and  requires  five  M-eoeff icients,  so  that  1 para- 
meters are  being  determined  from  lb  observations.  This  might  seem  likely  to  be  081 

disastrous,  and  it  is  not  surprising  that  the  sd  in  inclination  is  the  largest 


A 


I i 

among  the  85  orbits.  Comparison  with  the  more  accurate  and  more  reliable  orbit 
J shows  that,  although  the  values  of  i and  M,  differ  significantly,  the 
values  of  e,  , (u  + M^)  and  M,  are  consistent.  So  orbit  84,  though  poten- 
tially disastrous,  is  fairly  reliable.  it  • has  39  observations  and  covers 

2.5  days  before  decay:  six  M— coefficients  are  needed  and  M,.  is  exceptionally 

2 ^ 
large  (8.6  deg/day”).  All  analytical  orbit  determination  programs  break  down 

at  decay  because  the  perturbations  become  unlimited,  so  orbit  85  must  inevitably 

be  looked  on  with  suspicion.  Comparison  with  orbit  N shows  some  significant 

differences,  twice  the  sum  of  the  sd  on  inclination  and  five  times  on  eccentricity, 

but  good  agreement  on  . and  (u  + M() ) . So,  although  the  last  three  orbits  of 

Table  2 were  all  potentially  bad  (through  lack  of  observations  and  proximity  to 

decay),  they  emerge  remarkably  well  from  the  ust  of  comparison  with  the  (more 

accurate)  one-day  orbits. 

On  average  the  standard  deviations  on  orbits  83-85  of  Table  2 are  about 
three  times  larger  than  on  the  corresponding  one-day  orbits.  The  accuracy  of  the 
one-day  orbits  is  somewhat  better  than  those  of  the  main  orbits  without  Hewitt 
camera  observations,  though  not  so  good  as  the  orbits  with  Hewitt  camera 
observations . 

5 ANALYSIS  Of  VARIATION  IN  INC!. I!  . : ION  AND  ECCENTRICITY'  AT  RESONANCE 

If  we  accept  the  assumption  that  the  geopotential  can  be  expanded  in  a 
double  infinite  series  of  te  serai  harmonics,  'lumped'  harmonic  coefficients  of 
a particular  order  (linear  functions  of  individual  coefficients)  can  be  determined 
by  analysing  the  changes  that  occur  in  the  orbital  elements  of  satellites  vhi*.h 
experience  resonance  of  that  order.  The  satellite  1971—1 06A  was  appreciably 
perturbed  on  passing  through  14:1,  29:2  and  15:1  resonances,  and  the  effects  of 
these  resonances  on  the  orbital  inclination  and  eccentricity  have  been  evaluated. 

5 . 1 Theoretical  equations  for  3:a  resonance 

The  longitude-dependent  geopotential  at  an  exterior  point  (r,  P,  A)  may  be 
written  in  normalized  form'  as 


where  r is  the  distance  from  the  Earth's  centre,  9 is  co-latitude,  A is 
longitude  (positive  to  the  east),  is  the  gravitational  constant  for  the  Earth 

(398601  km  /s“) , R is  the  Earth’s  equatorial  radius  (6378,1  km),  pm(cos  8)  is 


12 


the  associated  Legendre  function  of  order  m and  degree  i , and  C„  and  S„ 

C,m  f.m 

are  the  normalized  tesseral  harmonic  coefficients.  The  normalizing  factor  N„ 

. . u 10 
is  given  by 

kt2  _ 2(2?,  + 1)U  - m)! 

N£m (F  + “ m)  ! * (3) 

The  rate  of  change  of  inclination  i caused  by  a relevant  pair  of  coef- 

. . - 1112 
ficients,  C„  and  S„  , near  6:a  resonance  may  be  written  ’ 

£m  Sim 


IF  = ~~^s~in  } ' (f)  hmPGSiPq(k  COS  1 ~ ~ jVexplj(Y^  " qw)f]  * 


- . . 12. 
where  Fn  is  Allan's  normalized  inclination  function  , G„  is  a function 
£mp  ’ £pq  ! ( 

of  eccentricity  e for  which  explicit  forms  have  been  derived  by  Gooding  , ■-# 

denotes  'real  part  of'  and  j = /-T  . The  resonance  angle  1>  is  defined  by  the 

equation 

$ = a(w  + M)  + 6(12  - v)  , (5) 


where  w is  the  argument  of  perigee,  M the  mean  anomaly,  12  the  right 

ascension  of  the  node  and  v the  sidereal  angle.  The  indices  y,  q,  k and  p 

in  equation  (3)  are  integers,  with  y taking  ihe  values  1,  2,  3 and  q 

the  values  0,  ±1,  ±2,  ; the  equations  linking  £,  m,  k and  p are*': 

ra  = yB;  k = ya  - q;  2p  = £ - k . 

At  B:a  resonance  the  m-suffix  of  a relevant  (C„  ,S„  ) pair  is  given 

Urn  Sim 

uniquely  by  the  choice  of  y . The  values  of  i to  be  taken  must  be  such  that 
i > m and  (S,  - k)  is  even.  The  successive  coefficients  which  arise  (for  given  y 
and  q ) may  usefully  be  gathered  together  in  a lumped  form  and  written  as'* 


V1  ci*k- 

i 


■ z 


where  i increases  in  steps  of  2 from  its  minimum  permissible  value  (2.  , and 

q k ...  U 

the  Q are  functions  of  inclination  that  can  be  taken  as  constant  for  a parti- 

q k 

cular  satellite,  and  = 1 . 

The  rate  of  change  of  eccentricity  e caused  by  the  (£,m)  harmonic  near 

0 u • 1 1 

p:a  resonance  can  be  written 
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<Ctm  - jSta)  e*p|j(Yt  - 1»)|J  . 


with  the  same  definitions  as  tor  equation  (4). 


(Ue) 


As  the  G t unctions''  are  ot  order  ^ \ , , , it  is  usually  found  that 

> in  ( i q | ) ! 

for  orbits  with  eccentricity  ^0.1  the  (t  ,q)  = (1,0)  terms  produce  the  most 
important  resonance  effects  on  the  inclination.  For  the  eccentricity,  the 

relative  importance  ot  ' he  terns  i-  largely  decided  by  the  value  of 

— G,  lu  - 1 (k  + q)e  which  is  of  order  lek  for  q = 0 , of  order  If  for 
e ipq 1 ’ , 

q = *1  and  of  order  J i e for  q - '2  . So  for  the  eccentricity  the  strongest 
effects  are  usually  caused  by  the  t>,q)  = (1.1)  and  (1,-1)  terms,  if  ie  < I . 

However,  these  rules  do  not  always  apply  and  1‘171-lObA  proves  to  be  some- 
what exceptional,  because  the  inclination  is  quite  close  to  the  critical  inclin- 
ation (63.4°)  and  the  variation  ot  uj  is  slow,  the  value  of  w being  about 
0 . S deg/day.  Consequently  the  (*1  - to)  , 'f  , and  (l1  + u>)  terms  in  equations  (4) 
and  (7)  are  difficult  to  separate;  when  all  three  are  included  in  the  fitting, 
the  correlations  between  them  tend  to  be  high,  and  there  is  a danger  that  the 
values  of  the  coefficients  obtained  will  have  large  standard  deviations.  So,  in 
analysing  the  three  resonances,  the  possibility  arises  of  dropping  one  of  the 
three  terms  to  improve  the  separation  of  the  coefficients.  This  had  to  be  tested 
on  all  three  resonances,  and  proved  a useful  stratagem  in  two  of  them. 

For  each  of  the  three  resonances  it  is  necessary  to  choose  the  time  interval 
over  which  the  analysis  is  to  be  made.  If  too  long  a time  interval  is  taken,  the 
orbital  parameter  fi  or  e)  will  not  be  appreciably  affected  by  the  resonance  near 
the  ends,  and  it  is  better  to  concentrate  the  analysis  in  the  region  where  the 
variations  are  strong.  The  danger  of  taking  too  short  an  interval  is  that  there 
will  not  be  enough  data  points.  A choice  between  these  extremes  must  be  made. 

In  practice  the  value  of  t is  used  as  a guide,  and  a range  of  values  of  i 
between  -b0  and  + hO  deg/day  is  regarded  as  the  outside  limit,  but  the  time  inter- 
val is  reduced  it  an  adequate  number  of  values  remain:  generally  it  is  advisable 
to  use  at  least  20  orbits. 

S . 2 I 4 1 h-order  (14:1)  resonance 

*3.2.1  Kquations  for  14:1  resonance 


The  most  important  terms  in  equation  (4)  for  14:1  resonance  are  those  with 
I , because  > = 2 terms  are  associated  with  harmonics  of  order  28  (m  = >|0, 


and  should  be  much  smaller  than  those  of  order  14.  Of  the  terms  with  > = I , 
those  with  q - 0,1  and  -I  are  likely  to  be  the  most  important,  since  terms 
with  q = +2  have  an  extra  e factor.  With  y = 1 , m = y 0 = 14  and 
k « >a  - q = I - q , and  concentrating  on  terms  with  (>  ,q)  = (1,0),  (1,1)  and 
(1,-1),  the  affixes  (q,k)  in  equations  (6)  are  (0,1),  (1,0)  and  (-1,2).  Writing 
only  tin-  three  terms  with  (>  ,q)  = (1,0),  (1,1)  and  ( 1 ,- 1 ) expl icitly  and  taking 

l I - e'  ) - I,  the  theoretical  variation  of  inclination  given  by  equation  (4)  may 

...  11-14 

be  written  tor  1 : I resonance  as 
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_q , k _q,k 

The  three  pairs  of  lumped  coefficients  C and  S appearing  in 

m m 

equation  (8)  may  be  written  in  terms  of  the  individual  geopotential  coefficients 

(C  ,S  ) as  indicated  in  equations  (b) . Explicitly,  with  the  Qq,k  expressed 

. Vm  - -q’k  14 

m terms  ol  the  K functions,  the  C are 
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.0,1 
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(II  ) 


and  similarly  for  S , on  replacing  C by  S throughout.  The  resonance  angle 
$ is  given  by  equation  (5)  with  a = I and  (3  = 14  . 
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The  most  important  terms  in  equation  (7),  which  gives  the  rate  of  change  of 
eccentricity,  are  those  with  (\  ,q)  = (1,1)  and  (1,-1),  but  for  consistency  with 
equation  (8)  the  (>,q)  = (1,0)  terms  are  also  given  explicitly.  The  theoretical 
equation  for  eccentricity  given  by  equation  (7)  may  therefore  be  written  for  14:1 


resonance  as 


- 15F. 


+ 1 1 F 


+ terms  n 


/ 0,1  0,1 

\ 

\ 1 4 Sin  * + C14 

cos  l1  j 

_i  ,o 

1 ,0 

' C,.  sin($  - io)  - 

S 1 4 C 

-1,2 

C,4  sin(f  + w) 

- S|  4 

<i«> IsUUI-'  , 

/ 1 _:t\-  i <9  ~ 

l(k  + 

|q  - i (k  + q)e2}  “ qu,)j|  ’ (1' 


_q,k  _q,k  il 

where  the  C and  S . are  given  by  equations  (9)  to  (II). 

I 4 I 4 

5.2.2  Analysis  of  inclination,  i 

Cosmos  4b2  passed  through  exact  14th-order  resonance  at  MJD  41659 
(1972  December  8).  The  effect  of  this  resonance  on  the  inclination  has  been 
analysed  over  a period  of  about  two  months  either  side  of  exact  resonance,  using 
the  THROE  computer  program  developed  by  Gooding'^’'',  which  fits  the  values  of  i 
with  equation  (4)  in  integrated  form.  During  this  time  there  were  27  values  of 
inclination  available  for  analysis,  10  values  from  the  PROP  orbits  in  Table  2 
and  17  values  from  orbits  supplied  by  the  US  Navy.  All  values  of  inclination 
were  cleared  of  lunisolar  and  zonal  harmonic  perturbations  using  the  computer 
program  PROD  with  one-day  integration  steps;  for  the  orbits  of  Table  2,  the 
J.j  9 perturbations  were  also  removed.  The  US  Navy  values  were  initially  given 
standard  deviations  of  0.003  and  the  PROP  values  were  given  their  quoted  standard 
deviations  from  Table  2.  Four  of  the  PROP  values  had  their  standard  deviations 
increased  to  0.0005°  to  allow  for  the  neglected  effect  of  Earth  tides;  one  of 
the  US  Navy  values  was  discarded  because  it  was  inaccurate  and  another,  that  at 
MJD  41612,  had  its  sd  increased  by  a factor  of  two. 

The  remaining  26  values  of  inclination  were  fitted  with  equation  (8)  using 

. . . . /-  -V1,k  • 

THROE,  ic  with  (> ,q)  = (1,0),  (1,1)  and  (1,-1).  The  coefficients  ^C,SJ|4  m 
equation  (8)  were  undetermined  in  this  fitting,  ie  the  standard  deviations  are 
of  the  same  order  as  the  values  of  the  coefficients.  This  was  to  he  expected 


089 


lb 


because  of  t lie  correlations  mentioned  earlier  (section  5 . I ) and  because  there 
were  only  2b  values  being  fitted,  scarcely  enough  to  allow  a good  determination 
of  the  seven  coefficients  (the  six  harmonic  coefficients  and  the  initial  value 
of  inclination). 

So  the  values  were  then  fitted  with  ((  ,q)  - (1,0)  only.  The  values  of 
lumped  harmonics  obtained  were 


9- 
10  0 


0, 


b2  ' 9 , 


9-°  * 1 

>0  s,,. 


* 19  , 


(13) 


with  > “ 1.6b  (t  is  as  defined  in  section  3).  A further  run,  with  (y,q)  = (1,0) 
and  (2,0),  offered  no  advantage:  the  (2,0)  coefficients  were  undetermined  and  t 
increased  slightly.  The  values  of  the  coefficients  in  equation  (13)  are  satis- 
factory in  that  they  provide  an  acceptable  fitting,  but  they  include  the  effects 

of  the  neglected  (1,1)  and  (1,-1)  terms.  Values  of  individual  lAth-order  eoef- 

0,1  0,1 


,14 


, . . '•  -9 

1971-IObA  should  be  approximately  9 \ 10  and  -25  x 10 


ficients  recent  lv  obtained’’  indicated  that  the  values  of  C.  and  S.,  for 

14  la 

These  values  are  of 

the  same  order  as  those  in  equation  (IS)  but  differ  sufficiently  to  suggest  that 

the  absorption  of  the  (1,1)  and  (1,-1)  terms  does  affect  the  values  obtained. 

9-° ’ 1 

Next  the  variation  ot  inclination  given  by  using  the  values  10  C.,  * 9 

' . 

and  10  S = -25  was  calculated  using  THROE.  It  was  found  that  the  variation 
1 4 

. . . . . o 

in  inclination  was  extremely  small,  never  more  than  *0.0008  . The  (1,0)  terms 
were  therefore  discarded  and  the  values  ol  inclination  fitted  with  ((  ,q)  = (1,1) 
and  (1,-1).  The  results  were  fairly  satisfactory,  with  c 1.5  and  standard 
deviations  about  \ of  the  values  of  the  coefficients. 

The  next  point  to  be  considered  was  the  possibility  that  (>,q)  = (1,2)  and 

/ - -V  1 

(1,-2)  terms  might  have  an  appreciable  effect.  Values  ol  \C,SJ  . and  yC.Sjj, 
were  calculated  using  the  values  of  individual  14th-order  coefficients  from 
Ref  14,  and  the  variation  of  inclination  due  to  the  (> ,q)  - (1,2)  and  (1,-2) 
terms  was  calculated  using  THROE.  These  terms  produced  an  appreciable  change  in 
inclination  at  resonance,  about  0.0015°.  The  raw  values  of  inclination  were  there 
fore  modified  by  subtracting  the  ettoct  ot  the  (1,2)  and  (1,-2)  terms  as  given  by 
this  calculation  and  a new  fitting  of  (>  ,q)  = (1,1)  and  (1,-1)  terms  was  made. 


OS" 


with  e = 1.48.  The  fitting  is  quite  good,  the  curve  being  shown  as  a full  line 
in  Pig  5. 

The  values  (14)  are  not  entirely  satisfactory  because  of  their  large 

standard  deviations,  which  stem  from  the  strong  correlations  between  coefficients  - 

_1 ,0  _-l,2 

for  example  the  correlation  between  C ^ and  is  -0.953.  So  a better 

solution  is  likely  to  be  possible  from  the  simultaneous  fit  of  inclination  and 
eccentricity  using  the  S1MRES  '*  computer  program.  This  will  be  discussed  after 
the  analysis  of  the  values  of  eccentricity. 

5.2.3  Analysis  of  eccentricity,  e 

The  effect  of  the  14:1  resonance  on  the  eccentricity  of  the  orbit  of 
1971  — 1 0f> A has  been  analysed  over  the  same  period  as  the  inclination,  ie  about  two 
months  either  side  of  exact  resonance,  using  the  same  27  orbits.  The  US  Navy 
values  were  given  standard  deviations  of  0.00008  and  the  PROP  values  were  given 
the  standard  deviations  quoted  in  Table  2.  Three  of  the  PROP  values  of  eccen- 
tricity had  their  standard  deviations  increased  to  0.000008  to  allow  for  the 
neglected  effect  of  Earth  tides,  and  the  US  Navy  value  dropped  from  the  inclin- 
ation analysis  because  of  inaccuracy  was  also  omitted  here,  for  the  same  reason. 

All  values  of  eccentricity  were  cleared  of  lunisolar  and  zonal  harmonic  pertur- 
bations using  the  PROD ' ^ computer  program  as  with  the  inclination  values. 

After  the  first  THROE  fitting,  with  equation  (12)  but  omitting  (\,q)  = (1,0) 
terms,  it  was  apparent  that  the  US  Navy  values  of  eccentricity  suffered  a bias 
relative  to  the  PROP  values.  This  mismatch  between  US  Navy  and  PROP  values  of  e 
has  arisen  previously'^,  and  was  attributed  to  an  inconsistency  in  restoring  the 
odd-zonal  harmonic  perturbation  to  the  US  Navy  values.  Since  the  perturbation  to 
e is  of  the  form  K sin  co  , where  K is  a constant,  the  simplest  procedure  for 
correcting  the  bias  is  to  subtract  an  expression  of  this  form  from  the  values  of 
e , the  value  of  K being  chosen  empirically  to  minimize  the  bias.  Here  the 
quantity  0.00015  sin  ui  was  subtracted  from  each  of  the  US  Navy  values  of  e . 
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1 

The  first  THROE  run  after  making  this  modification  showed  that  the  first 
PROP  value  of  e , at  MJD  41592,  fitted  badly  and  its  sd  was  increased  by  a 
factor  of  /TO.  The  resulting  THROE  run  for  (y,q)  = (1,1),  (1,-1)  gave  £ = 1.5 
and  values  of  the  coefficients  with  standard  deviations  not  much  smaller  than 
the  values.  Next  the  effects  of  the  terms  (y,q)  = (1,2),  (1,-2)  were  calculated 
using  values  of  the  individual  coefficients  from  Ref  14.  The  variation  in  e 
due  to  these  terms  was  appreciable,  about  0.00003,  so  the  values  of  e were 
cleared  of  this  perturbation,  and  refitted  with  THROE.  The  resulting  values  of 
the  coefficients  were: 


o_l,0 

10  C 1 4 

= - 139  + 

247  , 

oj.1 0 

10  S,4 

= 678  t 226 

q_-l.2 
10  C1 4 

= - 733  + 

154  , 

u.-1.2 

10  S,4 

> 

= -271  + 165 

with  £ = 1.50.  In  this  fitting  the  density  scale  height  H was  taken  as  42  km, 

the  value  appropriate  for  a height  ^H  above  perigee.  (Values  of  44  and  46  km  were 

also  tried,  but  were  less  satisfactory.)  The  fitting  is  quite  good,  the  curve 

being  shown  as  a full  line  in  Fig  6,  but  the  standard  deviations  of  the  values  are 

too  large  to  be  regarded  as  satisfactory,  again  because  of  the  high  correlation 

.1,0  _H,2 

between  coefficients  - for  example  the  correlation  between  S,,  and  C,,  is 

14  14 

0.955.  Therefore  a simultaneous  fitting  with  i is  required  and  is  discussed  in 
the  following  section. 

5.2.4  Inclination  and  eccentricity  fitted  simultaneously 

The  values  of  inclination  and  eccentricity  fitted  separately  by  THROE  were 
next  fitted  simultaneously  using  the  computer  program  SIMRES  ' ' . This  program 
combines  the  results  from  a number  of  THROE  runs  and  produces  a single  set  of 
coefficients  to  fit  the  data.  The  program  allows  a choice  of  weighting,  so  that 
the  contributing  THROE  runs  can  be  given  more  or  less  weight  according  to  their 
accuracy  of  fit,  indicated  by  the  value  of  e . 

Here  the  final  THROE  runs  for  i and  e were  combined  and  no  weighting 
factor  was  applied  because  e was  1.5  in  both  contributing  THROE  runs.  The 
SIMRES  fittings  are  shown  in  Figs  5 and  6 by  broken  lines  and  the  values  of  the 
coefficients  given  by  SIMRES  are: 
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10  S, , =-137'  35 
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10  S,  , = b6  * 42 
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The  standard  deviations  ot"  the  values  in  (lb)  are  considerably  smaller  than 
in  (14)  or  (15),  and  the  values  should  also  be  more  reliable  because  the  number 
of  values  being  fitted  is  doubled.  This  conclusion  is  confirmed  by  Figs  5 and  6, 
which  show  that  for  both  i and  e the  SIMRES  fitting  looks  just  as  good  as  - 
or  perhaps  better  than  - the  individual  fittings.  The  values  (lb)  differ  from 
the  corresponding  values  in  (14)  and  (15)  by  between  0.5  and  3.2  times  the  sum 
of  the  standard  deviations.  The  individual  differences  for  each  coefficient, 
expressed  as  multiples  of  the  sum  of  the  standard  deviations,  are  as  follows: 
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S14 

L 14 

S14 

1 .5 

2.8 

2.5 

1 .8 

0.5 

3.1 

3.2 

1 .b 

1,0  -1,2 

From  this  table  it  might  be  expected  that  the  values  of  and  S , ^ 

be  more  reliable  than  the  other  two. 


wou  1 d 


In  the  recent  evaluation  of  individual  14th-order  harmonies  in  the  geo- 
14  ... 

potential  , the  lumped  harmonic  values  in  equation  (.lb)  were  used  in  the  solutions 

1,0  _-l,2 

for  harmonics  of  order  14  and  even  degree.  It  was  found  that  the  C.,  and  S,, 

14  14 

values  in  equation  (lb)  made  a useful  contribution  to  the  solution,  yielding 

weighted  residuals  of  0.04  and  0.71  respectively  in  the  five-coefficient  solution 

1,0  -1,2 

quoted  in  Table  7 of  Ref  14.  The  S, , and  C,.  values  in  equation  (lb),  how- 

14  I 4 

ever,  did  not  fit  so  well  and  their  standard  deviations  were  increased  by  factors 
of  5 and  10  respectively:  even  then,  the  weighted  residuals  for  these  two  coef- 
ficients in  the  five-coefficient  solution  quoted  in  Table  7 of  Ref  14  are  -1.17 

_1,0 

and  -1.14  respectively.  This  confirms  the  indication  that  the  values  of  0,, 

1,2  14 

and  S.,  in  equation  (lb)  are  more  reliable  than  the  other  two. 

I 4 

5.3  24:2  resonance 

Tin’  changes  in  inclination  and  eccentricity  at  29:2  resonance  are  expected 
to  be  only  k as  large  as  at  14:1  resonance,  because  the  values  of  the 


are 


20 


likely  to  be  only  ^ as  large,  and  the  resonance  is  faster  because  $ — aM  is 
twice  as  large.  So  the  chances  of  successful  analysis  are  poorer.  However,  the 
analysis  needs  to  be  made  to  assess  the  overall  change  in  i at  resonance. 

5.3.1  Equations  for  29:2  resonance 

The  most  important  terms  in  equation  (4)  for  29:2  resonance  are  those  with 
y = 1 and  q = 0,  1 and  -I.  The  y = 2 terms  are  associated  with  harmonics  of  order 
58  (m  = yg),  which  should  be  much  smaller  than  those  of  order  29;  and  terms  with 
q = ±2  have  an  extra  e factor. 

For  the  29:2  resonance  with  y = I , m = > f?  = 29  and  k = ya-q  = 2-  q,  so 

that  the  affixes  (q,k)  in  equation  (6)  are  (0,2),  (1,1)  and  (-1,3)  when  q = 0, 

I and  -1  resp  ctively.  Writing  only  the  three  terms  with  (> ,q)  = (1,0),  (1,1) 

and  (1,-1)  explicitly,  the  theoretical  variation  of  inclination  given  by 

13  18 

equation  (4)  may  be  written  for  29:2  resonance  ’ as 
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The  three  pairs  of  lumped  coefficients  C and  S appearing  in 

mm 

equation  ( 1 7)  may  be  written  in  terms  of  the  individual  geopotential  coefficients 


(Sm,Sem) 


• ok 

as  indicated  in  equations  (6).  Explicitly,  with  the  Q,  expressed 


_q,k 


in  terms  of  the  F functions,  the  C are 
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and  similarly  for  S , on  replacing  C by  S throughout.  The  resonance  angle 
T is  given  by  equation  (5)  with  a = 2 and  t>  - 29. 

For  the  29:2  resonance,  the  theoretical  variation  of  eccentricity  given  by 

1 8 

equation  (7)  may  be  written  in  terms 
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Three  terms  are  given  explicitly  in  equation  (21),  those  with  (y,q)  = (1,0), 

(1,1)  and  (1,-1).  The  main  terras  are  expected  to  be  those  with  (y,q)  = (1,1) 
and  (1,-1)  but  the  (y,q)  = (1,0)  term  is  also  given  explicitly  for  consistency 
with  equation  (17). 

3.3.2  Analysis  of  inclination,  i 

The  inclination  of  Cosmos  462  was  analysed  at  the  time  of  29:2  resonance, 
using  21  values  of  inclination,  seven  values  being  from  the  PROT  orbits  in 
Table  2 and  14  values  from  orbits  supplied  by  the  US  Navy.  The  analysis  extended 
to  nearly  two  months  either  side  of  exact  29:2  resonance,  which  occurred  at 
MJP  42012  (1973  November  2t>).  All  values  of  inclination  were  cleared  of  1>  ni- 
solar  and  zonal  harmonic  perturbations  using  the  PROD  * ^ program  with  one-day 
integration  steps  and,  for  the  seven  values  from  Table  2,  the  perturbation. 

The  US  Navy  values  were  given  standard  deviations  of  0.003°  and  the  PROP  values 
were  given  their  quoted  standard  deviations  in  Table  2. 

The  21  values  of  inclination  were  then  fitted  with  equation  (17)  using 

. . I-  -\q,k 

THROE,  ic  with  (>,q)  = (1,0),  (1,1)  and  (1,-1).  All  the  (c,Sj)9  coefficients 


in  equation  (.1 7 1 were  undetermined  in  this  fitting,  but  t lie  tit  did  reveal  that 

the  US  Navy  values  were  all  too  large:  discrepancies  of  this  kind  have  been 

17, 14,20  . ... 

found  previously  and  are  due  to  differences  in  definition  in  the  l'S  Navy 

orbits.  The  difficulty  is  avoided  by  making  a (constant)  empirical  change  in 

the  VS  Navy  values  and  here  they  were  all  reduced  by  0.004  . The  standard 

deviations  on  one  PROP  value  of  itn  lination,  at  MJU  42018,  and  one  VS  Navy,  at 

MJl)  42060,  were  increased  by  a factor  of  two.  The  program  was  then  re-run,  but 

with  similar  results  - all  seven  coefficients  were  undetermined,  as  was  expected 

because  of  the  correlations  mentioned  earlier  (see  section  8.1)  and  also  because 

there  were  only  21  values  being  fitted. 

The  values  were  then  fitted  with  (>,q)  ■*  (1,0)  only,  because  at  this 
resonance  e = 0.065  and  therefore  the  (1,1)  and  (1,-1)  terms  should  not  have  so 
much  effect.  The  values  obtained  were 

0 2 o ’ 

* - o_  » 

10  CH)  ■=  -90  + 74  . 10  S = -127*34  , (22) 


with  i = 0.689.  These  values  were  obtained  after  the  M,  values  on  the  orbits 
had  been  changed  (see  section  5.3.1),  The  fitted  curve  is  shown  as  a full  line 
in  Kig  7.  The  values  in  equation  (22)  should  bo  useful  in  the  future  in  determin- 
ing the  individual  coefficients  of  2°th  order  and  even  degree. 

Even  though  the  numerical  values  ol  the  lumped  harmonics  were  not  particu- 
larly accurate,  the  fitting  is  useful  in  showing  t li.it  there  was  an  overall  de- 

o ... 

crease  ot  0.002  in  passing  through  2l|:2  resonance  and  m providing  end  points  tor 

the  analysis  ot  atmospheric  rotation  (see  section  6). 


5.3.1  Analysis  of  eccentricity,  e 
* - 


The  effect  of  the  29:2  resonance  on  the  eccentricity 
values  of  e from  the  same  21  orbits.  All  values  of  e 
solar  and  zonal  harmonic  perturbations  using  PROo'*  . 1'he 
the  standard  deviations  quoted  in  Table  2 and  the  VS  Navy 
standard  deviations  ot  0.00008.  l’he  density  scale  height 

After  the  lirst  THROE  fitting  with  equation  (211, 
(1,1)  and  (1,-1),  it  was  again  apparent  that  the  l'S  Naw  v 
a bias  relative  to  the  1’KOI’  values.  1'his  bias  was  correct 
section  5.2.3)  by  subtracting  an  expression  ol  the  form  k 
a constant  chosen  empirically  tv’  minimize 
The  subsequent  THROE  fitting  witli  (> ,ql 
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undetermined  values  of  the  coefficients  and  they  also  had  high  correlations:  it 
was  obviously  not  possible  to  determine  seven  coefficients  from  21  values. 

In  the  course  of  the  fitting  of  e , a source  of  error  in  the  technique  for 
using  the  THROE  program  became  apparent.  The  correction  of  eccentricity  to  allow 
for  the  effects  of  atmospheric  drag  between  one  orbit  and  the  next  is  calculated 
in  THROE  using  the  value  of  M,  for  the  first  of  the  two  orbits,  which  is 
assumed  to  apply  over  the  time  interval  between  them.  However,  it  is  not  correct 
to  assume  that  the  value  of  M,,  for  the  first  orbit  remains  valid  over  the  time 
interval  between  the  two  orbits,  often  seven  days:  for  the  PROP  orbits  the  value 
of  M,  usually  applies  over  about  three  days  either  side  of  epoch,  and  for  the 
US  Navy  orbits  for  approximately  five  days  before  epoch.  So  the  technique  hither- 
to used  was  modified,  and  the  value  of  M„  at  the  nth  epoch  was  taken  as 

<(M,)  . - (M.)  >/2(t  , - t I . This  ensures  that  the  integrated  effect  of 

( I n+ 1 1 n )/  \ n+ 1 n / 

air  drag  between  epoch  t and  epoch  t j is  correctly  represented.  This 
correction  had  an  important  effect  in  the  eccentricity  runs,  but  no  significant 
effect  for  inclination,  where  the  atmospheric  corrections  are  very  small;  however, 
the  values  of  M,,  were  changed  for  i , for  consistency  between  the  fittings. 

After  changing  the  values  of  M0  , the  values  of  e to  be  fitted,  when 
adjusted  for  the  effects  of  atmospheric  drag  by  THROE,  showed  very  little  varia- 
tion, all  the  values  being  within  0.00015  of  the  value  0.06900.  So  the  chances 
of  obtaining  good  lumped  values  of  29th-order  harmonics  were  unpromising.  Fitting 
with  (>  ,q)  = (.1,1)  and  (1,-1)  gave  quite  undetermined  values,  so  it  was  necessary 
to  fit  with  (>  ,q)  = (1,1)  only,  or  (> ,q)  = (1,-1)  only.  These  alternative  fit- 
tings were  very  similar  with  c = 1.33  and  1.30  respectively,  but  the  latter  seems 
marginally  preferable  because  it  fits  the  PROP  values  of  eccentricity  better. 

The  (y,q)  = (1,-1)  fitting  is  shown  in  Fig  8. 

The  values  obtained  for  the  coefficients  in  this  (1,-1)  fitting  are 

— 9 — o . — 1,3 

(907  * 325)  ' 10  and  (-463  ’163)  x 10  : these  are  nominally  values  of  C,1c) 

_-l ,3  _ 1 , 1 _ 1 , 1 

and  respectively,  but  they  include  the  effects  of  C,lt)  and  S 

J,  1 - J,1 

However  the  results,  even  if  the  interference  from  and  is  slight, 

are  not  accurate  enough  to  give  numerical  values  useful  in  determining  individual 
29th-order  coefficients.  In  other  words,  the  variation  in  eccentricity  over  the 
region  of  29:2  resonance  is  too  small  to  yield  good  values  of  lumped  coefficients 
because  the  orbit  of  1971—1 06A  is  strongly  affected  by  drag  and  the  effect  of 
089  resonance  does  not  last  long  enough  for  an  appreciable  change  in  e to  build  up. 
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Since  the  change  in  e is  not  large  enough  to  give  values  of  lumped 
coefficients,  there  was  no  purpose  in  making  a SIMRES  fitting  of  i and  e 
together. 

5.4  15th-order  (15:1)  resonance 

5.4.1  Equations  for  15:1  resonance 

As  for  the  two  preceding  resonances,  the  most  important  terms  in  equation 
(4)  for  15:1  resonance  are  likely  to  be  those  with  y = I,  because  y = 2 terms 
are  associated  with  harmonics  of  order  30  (m  = yb) , and  should  be  considerably 
smaller  than  those  of  order  15.  Again  the  most  important  terms  with  y = 1 are 
those  with  q = 0,1  and  -1,  since  q = 2 terms  have  an  extra  e factor.  Therefore, 
with  y = 1,  m = 15  and  k = ya  - q = 1 - q,  and  concentrating  on  terms  with 

(y,q)  = (1.0),  (1,1)  and  (1,-1),  the  affixes  (q,k)  in  equations  (6)  are  (0,1), 
(1,0)  and  (-1,2).  So,  writing  only  the  three  terms  with  (y,q)  = (1,0),  (1,1) 
and  (1,-1)  explicitly,  equation  (4)  giving  the  theoretical  variation  of  inclin- 


ation may  be  written  for  15:1  resonance 
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The  tnree  pairs  of  lumped  coefficients  and  appearing  in 

equation  (23)  may  be  written  in  terms  of  the  individual  geopotential  coefficients 


as  indicated  in  equation  (6).  Explicitly,  with  the  Q,  expressed 
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in  terms  of  the  F functions,  the  C are 
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and  similarly  for  S , on  replacing  C by  S throughout.  The  resonance  angle 
i is  given  by  equation  (5)  with  a = 1 and  8=15. 

For  the  15:1  resonance,  the  theoretical  variation  of  eccentricity  given  by 

/q,k  _q,k\  j 2 

— /-.N  — . * .u e (c  s ) — ’ 

\ m m / 


equation  (7)  may  be  written  in  terms  of  the  same 
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Three  terms  are  given  explicitly  in  equation  (27),  those  with  (y,q)  = (1,0),  (1,1) 
and  (1,-1).  The  main  terms,  as  with  the  previous  two  resonances,  are  expected  to 
be  those  with  (y,q)  = (1,1)  and  (1,-1),  but  the  (y,q)  = (1,0)  term  is  given  for 
consistency  with  equation  (23). 

5.4.2  Analysis  of  inclination,  i 

At  the  time  of  )5th-order  resonance,  nine  PROP  orbits  from  Table  2 and  13 
US  Navy  orbits  were  available  for  analysis.  The  analysis  covers  a period  of 
approximately  two  months  before  and  one  month  after  the  date  of  exact  15:1  reson- 
ance, MJD  42302  (1974  September  12).  All  values  of  inclination  were  cleared  of 
lunisolar  and  zonal  harmonic  perturbations  using  the  PROd'^  program  with  one-day 
integration  steps;  and  the  nine  values  from  Table  2 were  cleared  of  the 
perturbation.  The  US  Navy  values  were  given  standard  deviations  of  0.003°  and 
the  PROP  values  were  given  their  quoted  standard  deviations  in  Table  2. 
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The  22  values  of  inclination  were  fitted  with  equation  (23)  using  THROE, 

q,k 

ie  with  (y ,q)  = (1,0),  (1,1)  and  (1,-1).  The  coefficients  (C,S)j^  in 
equation  (23)  were  undetermined  in  this  fitting,  so  the  values  were  then  fitted 
with  just  (y,q)  = (1,0).  For  this  fitting  the  US  Navy  value  at  MJD  42333  was 
dropped  because  it  was  inaccurate  and  the  values  on  the  orbits  were  changed 

as  for  the  29:2  resonance  (see  section  5.3.3).  The  values  obtained  were 

o_0.>  o-O.' 

10  C[5  = -36+21  , 10  S15  = 9 ± 17  , (28) 


with  e = 1.102.  The  fitted  curve  is  shown  as  a full  line  in  Fig  9.  These  values 
of  (C,S)°!‘  are  quite  small  and  consequently  not  well  determined.  Previous 

1 J 2 ] Q 

evaluations  of  1 5th-order  coefficients  indicate  that  at  inclinations  near  65 

_0,1  _0,1  _g 

the  values  of  both  C j ^ and  S|^  are  near  zero,  and  of  order  (0  ± 20)  x 10  ; 

-0,1  _0,1  0 . 

the  variation  of  and  near  65  is  not  very  well  defined  in  Ref  21 

because  the  only  value  available  near  65°  - Cosmos  387  at  62.9°  - may  not  be 

entirely  reliable,  as  the  data  on  which  it  is  based  are  very  limited. 

As  the  (C,S)°’]  coefficients  are  small  for  J971-106A,  it  was  worth  fitting 
the  inclination  values  with  (y,q)  = (1,1)  and  (1,-1).  The  values  obtained  were 


9-1  0 

oJ 

196  ± 

- 

,0  C15  = 

q_->  *2 

29  ± 83  , 

10  S15  = 

q_~»  .2 

142 

> 

10  C, 5 = 

116  ± 124  , 

10  S 1 5 = 

145  ± 

97 

a 

(29) 


with  e = 1 .096.  These  results  were  then  used  with  the  results  for  eccentricity 
in  a SIMRES  fitting  (section  5.4.4). 

5.4.3  Analysis  of  eccentricity,  e 

The  same  21  orbits  used  in  analysing  the  inclination  were  analysed  to 
determine  the  effect  of  15:1  resonance  on  the  eccentricity:  the  US  Navy  value 
at  MJD  42333  was  dropped  again  because  it  was  inaccurate.  All  the  21  values  of 
e were  cleared  of  lunisolar  and  zonal  harmonic  perturbations  using  PROd'^.  The 
PROP  values  were  given  the  standard  deviations  quoted  in  Table  2 and  the  US  Navy 
values  were  given  standard  deviations  of  0.00008.  The  density  scale  height  was 
taken  as  35.5  km. 
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At' tot  the  tirst  TIIROT  titling  w i tit  equation  (27)  it  was  again  apparent  tliat 
the  US  Navy  values  suffered  a bias  relative  to  t lie  l'KOP  values  (see  seel  ion  5.7.  1). 
I’o  overcome  this  bias,  the  quantity  0.000.’  sin  10  was  subtracted  1 rom  each  US  Navy 

value  ot  eeeent r i e i t y . In  order  to  obtain  the  right  adjustment  lor  the  elteets 

ot  atmospheric  drag  by  THROK , the  values  ot  M,  were  changed  using  the  technique 
described  in  sect  ion  5.  1.  t. 

Alter  the  tirst  tilting  with  THROE,  with  (>  ,q)  - (1,0),  (1,1)  and  (1,-1), 

it  was  apparent  that  seven  cocfticients  could  not  hi'  determined  from  7 1 values  ot 

eccentricity.  As  the  (1,0)  terms  have  very  little  effect  on  eccentricity,  the 
values  were  next  fitted  with  (>  ,q)  (1,1)  and  (1,-1).  The  values  of  the  lumped 

coefficients  obtained  were 


i 24  • bA  , 


q- 1 o 

10  S ’ - 28  * 1 04 


- 1 54  • 8l>  , 10  S,,  - - b ' 54 
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with  t - I . 4 b h . The  fitted  curve  is  shown  as  a full  line  in  Tig  10.  The  values 
of  all  the  coefficients  in  equation  (TO)  agree  with  those  from  equation  (2l>)  to 
within  1.  1 times  the  stun  ot  their  standard  deviations.  So  the  prospect  of 
achieving  a better  fit  with  S1MKKS  is  good. 

5.4.4  Inclination  and  eeeent  r i e i t y lilted  simultaneously 

The  values  of  inclination  and  eccentricity  fitted  separately  by  THROK  with 
(1  ,q)  ” (1,1)  and  (1,-1)  were  next  fitted  simultaneously  using  the  computer 
program  SIMRKS* * . As  explained  in  section  5.2.4,  this  program  allows  a choice 
ot  weighting.  In  the  tirst  fitting,  i and  e were  given  equal  weight  and  in 
the  second  e was  degraded  by  a factor  equal  to  the  ratio  of  the  final  values  ot 
t on  the  THROK  fittings,  namely  1.140  (»  I .4b9/ 1 .Odb) . The  second  fitting  is 
the  more  logical,  and  gave  lower  standard  deviations  for  all  four  lumped  coef- 
ficients; so  it  was  preferred.  The  tit  tings  are  shown  in  Tigs  l>  and  10  by 
broken  lines,  and  it  is  very  satisfactory  that  these  are  just  as  good  as  the 
separate  fittings  (unbroken  lines).  The  values  of  t ho  lumped  coefficients  given 
bv  S I MKTS  are 
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It  can  be  seen  ^hat  the  standard  deviations  here  are  much  lower  than  those  in 
equations  (29)  and  (30),  so  the  combined  fitting  is  well  worthwhile.  The  values 
in  equations  (31)  differ  from  the  corresponding  values  in  equations  (29)  by  less 
than  1.7  times  the  sum  of  the  standard  deviations  and  from  the  corresponding 
values  in  equations  (30)  by  less  than  0.9  times  the  sum  of  the  standard  deviations. 
The  individual  differences  for  each  coefficient  expressed  as  multiples  of  the  sum 
of  the  standard  deviations  are  as  follows: 
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-1  .2 
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The  lumped  coefficients  in  equations  (31)  should  be  useful  in  a future 

determination  of  the  individual  coefficients  of  the  15th-order  and  even  degree. 

22 

In  a previous  determination  “ there  were  no  reliable  values  for  inclinations 
near  65  . 

6 ATMOSPHERIC  ROTATION 

6 . I Introduction 

The  upper  atmosphere  is  rotating  at  approximately  the  same  rate  as  the 
Earth:  therefore  the  aerodynamic  force  acting  on  a satellite  has  a component 

perpendicular  to  the  orbit,  and  this  has  the  effect  of  reducing  the  inclination, 
i , of  the  orbit  during  the  course  of  the  satellite's  life.  For  a high-drag 
satellite,  atmospheric  rotation  is  the  most  important  force  perturbing  i ; so 
if  the  change  in  i is  accurately  measured,  and  other  perturbations  are  removed, 
the  rotation  rate  of  the  upper  atmosphere  in  the  region  near  the  satellite's 
perigee  can  be  determined,  ie  the  zonal  (west-to-east)  wind  speed  near  perigee 
can  be  evaluated. 

The  perturbations  to  be  removed  are  caused  by  (1)  lunisolar  gravitational 
attraction,  (2)  zonal  harmonics  in  the  geopotential , (3)  the  J.v,  tesseral 
harmonic,  (4)  any  change  in  i due  to  resonance  and  (5)  the  change  in  inclina- 
tion due  to  meridional  winds.  The  values  of  inclination  are  cleared  of  the 

I 6 

effects  of  (1)  and  (2)  using  the  computer  program  PROD  with  one-day  integration 
steps.  The  effect  of  perturbation  (3)  is  removed  by  calculation  of  its  numerical 
value  for  each  PROP  orbit.  The  resonances  have  already  been  analysed  in 
section  5,  and  the  effects  of  meridional  winds  and  other  perturbations  are 
discussed  in  section  6.3. 
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6. 2 Theory 

The  change  Ai  in  the  inclination  of  a satellite's  orbit,  due  to  atmospheric 

rotation,  is  given  in  terms  of  the  change  AT  in  the  satellite's  orbital  period, 

23  d 

due  to  drag,  by  the  equation 
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(32) 


In  equation  (32),  which  applies  for  an  eccentricity  of  0.2  or  less,  A is  the 
angular  velocity  of  the  atmosphere  (about  the  Earth's  axis)  in  the  region  near 
perigee,  divided  by  the  Earth's  angular  velocity:  thus  A is  strictly  non- 
dimensional,  but  can  conveniently  be  expressed  in  rev/day  because  the  Earth's 
rotation  rate  is  1.0  rev/day.  In  equation  (32)  fhe  change  in  orbital  period  AT^ 
is  expressed  as  a fraction  of  a day;  the  term  in  cos  4oi  , which  was  neglected 
in  Ref  23,  has  been  restored. 

Other  parameters  in  equation  (32)  are  as  follows:  the  I are  the  Bessel 

n 

functions  of  the  first  kind  and  imaginary  argument,  of  degree  n and  argument 
z = ae/H  , where  H is  the  atmospheric  density  scale  height.  The  parameter  c 
takes  account  of  atmospheric  oblateness  and  is  given  by  c = |e'a(l  - e)  sin^i|/2H  , 
where  e'  is  the  ellipticity  of  the  atmosphere,  taken  the  same  as  that  of  the 
Earth,  0.00335.  The  factor  F is  given  by  /F  = jl  - a(l  - e)w  cos  i|/V^  , 
where  is  the  satellite's  velocity  at  perigee  and  w is  the  angular  velocity 

of  the  atmosphere  near  perigee.  Usually  /f"  has  a value  between  0.95  and  1.05. 

For  1971-106A,  at  inclination  65.7°  with  perigee  height  near  200  km, 
v/F  — 1 - 0.0025A  . 
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The  change  in  inclination  due  to  meridional  winds  is  ttiven  by 
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2.2. 

where  K = sin~i/(l  + cos'i)  and  m is  the  south-to-north  atmospheric  rotation 

rate,  that  is  the  south-to-north  wind  speed  (m/s)  divided  by  rw  where  r is 

the  distance  from  Earth's  centre  (m)  and  w is  the  angular  velocity  of  the 

Earth  (72.7  * 10  rad/s);  so  rw  = 6580000  * 72.7  * 10  = 478  m/ s if  perigee 

h. 

height  is  near  200  km. 

6 . 3 Fitting  of  theoretical  curves  to  the  inclination 

The  values  of  orbital  inclination  from  the  85  orbits  of  Table  2 and  the 
15  orbits  of  Table  4 have  been  cleared  of  the  effects  of  lunisolar,  zonal  harmonic 
and  J0 ^ tesseral  harmonic  perturbations,  and  the  values  have  been  plotted  as 
circles  in  Figs  II  and  12.  The  15  values  of  Fig  12  can  be  regarded  as  super- 
seding the  last  three  values  of  Fig  11. 

Perturbations  in  inclination  due  to  Earth  and  ocean  tides  and  solar  radia- 
tion pressure  have  not  been  taken  into  account.  The  tidal  perturbation  in  i 

2C 

for  Geos  1,  a satellite  with  an  orbital  inclination  similar  to  1971-I06A,  amounts'" 

to  ±0.0004°.  To  allow  for  the  neglect  of  this  effect,  all  the  standard  deviations 

less  than  0.0005°  have  been  increased  to  0.0005°.  The  effect  of  solar  radiation 

pressure  is  very  small:  its  effect  on  Explorer  24,  a satellite  in  a comparable 

26 

orbit  but  of  much  higher  area-to-mass  ratio  (which  was  analysed  by  Slowey  ), 
shows  that  the  effect  of  solar  radiation  pressure  on  i for  Cosmos  462  is  likely 
to  be  of  order  0.00001°. 

The  theoretical  variation  of  i due  to  atmospheric  rotation  and  meridional 
winds  can  be  determined  for  a series  of  values  of  A and  p , using  a computer 
program  (ROTATM)  based  on  equations  (32)  and  (33).  The  theoretical  variation  of 
i was  calculated  for  values  of  A between  0.8  and  1.4  at  intervals  of  0.1,  with 
P = 0 and  0.2,  every  20  days  throughout  the  satellite's  life  until  the  last  15 
days,  when  the  variation  of  i was  evaluated  daily.  The  values  of  i in  Figs  11 
and  12  were  then  fitted  with  the  best  theoretical  curves.  Fig  II  divides  natur- 
ally into  four  sections  separated  by  the  perturbations  in  inclination  at  14:1, 

29:2  and  15:1  resonances.  These  resonances  have  been  analysed  in  section  5 and 
the  change  in  inclination  due  to  resonance  is  available  for  each  resonance.  At 
the  end  of  each  resonance  period  the  curve  for  atmospheric  rotation  is  started  at 
the  value  of  i given  at  the  end  of  resonance  in  the  resonance  analysis  (see 
Figs  5,  7 and  9),  with  atmospheric  rotation  perturbations  restored.  There  is  a 
fifth  section  to  be  fitted  with  a theoretical  curve,  for  the  last  15  days  when 
the  daily  orbits  are  available  (Fig  12). 
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I'lic  variations  in  i due  to  meridional  rotation,  south-to-north  winds, 

ire  proportional  to  cos  m (see  equation  (33))  and  therefore  the  effect  usually 

cancels  out  over  one  cycle  of  u>  for  a fixed  value  of  p . However,  the  value 

27-29  . 29-32 

ol  itscli  varies  and  both  theoretical  studies  and  experimental  results 

suggest  that  has  an  approximately  sinusoidal  variation  during  the  course  of 

23  hours;  with  maximum  wind  towards  the  equator  at  about  02  h local  time  and 
maximum  wind  away  from  equator  at  about  14  h.  For  1971—1 0fc>A  the  variation  in  w 
is  slow  (one  cycle  in  1.8  yr  approximately,  see  Fig  3)  and  the  local-time  varia- 
tions, indicated  at  the  top  of  Fig  11,  are  much  faster,  one  cycle  every  90  days 
approximately.  The  first  four  sections  in  Fig  11  are  all  averaged  in  local  time, 
covering  about  200  days.  The  effect  of  p will  therefore  tend  to  cancel  out, 
and  the  fittings  were  made  with  the  p = 0 curves.  It  is  possible  that  the  mean 
value  of  p is  non-zero,  but  the  effect  is  generally  small:  on  the  first  section, 
for  example,  the  effect  of  p = 0.05  throughout  would  be  to  change  i by  0.0005°. 

So  the  values  of  inclination  are  fitted  by  choosing  the  best  values  of  A 
assuming  p = 0 on  the  first  four  sections.  After  allowance  has  been  made  for  the 
breaks  at  resonance,  the  best  fit  between  theory  and  observations  for  the  four 
sections  of  Fig  II  was  obtained  with  values  of  atmospheric  rotation  A of  1 . 1 , 

0.9,  1.0  and  1.0  respectively. 

In  order  to  assess  the  likely  errors  in  these  values  of  A , a realistic 
estimate  of  the  likely  errors  in  i at  the  beginning  (i„)  and  end  (i.J  of  each 

J 2 2 / B E 

section,  a and  o , was  made,  and  then  Arc  + a /(i  ~ i„)  was  taken  as 

B fci  B E / B E 

the  standard  deviation  in  A . For  the  first  three  sections  of  the  curve  in  Fig  11 
the  errors  o and  are  estimated  as  being  between  0.0005C  and  0.0007°,  and 

this  gives  errors  in  A of  near  0.05  on  all  three  sections.  For  the  fourth 
section  the  error  in  A is  0.02  on  the  same  basis.  This  value  is  so  small  that 
other  errors,  hitherto  neglected,  may  become  significant:  there  is  a possible 
error  of  about  0.01  from  the  neglect  of  O(c^)  terms  in  equation  (32),  and  possibly 
errors  approaching  0.01  from  errors  in  the  assumed  atmospheric  ellipiticity . 

There  may  also  be  a small  error,  not  more  than  0.02,  from  neglect  of  the  effects 
of  31:2  resonance. 

There  is  a possibility  that  the  atmospheric  winds  depend  on  season,  so  the 
latitude  and  'seasonal  bias'  of  perigee  are  indicated  below  the  curve  in  Fig  II. 

In  fitting  the  values  of  inclination  in  the  last  15  days  of  the  life.  Fig  12, 
the  effect  of  meridional  winds  should  be  included.  The  season  is  near  equinox 

and  the  perigee  latitude  is  54-b2  N:  in  these  conditions,  the  prevailing  winds 
27-32 

are  “ from  north  to  south  of  order  100  m/ s at  local  times  from  23  h through 
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midnight  to  05  h,  with  only  very  weak  winds  from  17-23  h local  time.  So  ii  is 
appropriate  in  Fig  12  to  take  g = -0.2  up  to  March  30,  and  n = 0 thereafter. 

With  these  values  of  p , the  best  fit  is  obtained  with  A = 1.1:  the  curve  fits 
all  the  15  values  in  Fig  12  to  within  1.5  sd,  and,  taking  o and  Oj.  as  slightly 
less  than  0.001  , the  sd  in  A is  assessed  as  0.1. 

b . 4 Results 

33 

Previous  studies  of  upper-atmosphere  zonal  winds  have  shown  that  the  value 
of  A varies  with  both  height  and  local  time.  Fig  13  of  Ref  33  gives  three  curves 
for  the  variation  of  A with  height,  for  evening  (18-24  h),  morning  (04-12  h) 
and  average  values  of  local  time.  A revised  version  of  this  diagram  is  given  as 
Fig  2 of  Ref  34. 

Tab 1 e 5 

Values  of  atmospheric  rotation  rate,  A , obtained  from  1971-1 0b A 


Date 

Height 

km 

Local 

time 

1 Seasonal 
bias’ 

21  Jan  - 2 Oct  1972 

253 

Av 

Winter 

1.1  0.05 

30  Jan  - 30  Sep  1973 

250 

A v 

Summe  r 

0.9  i 0.05 

20  Jan  - 14  Jul  1974 

239 

Av 

Average 

1.0  i 0.05 

24  Oct  1974  - 19  Mar  1975 

218 

Av 

Average 

1 .0  1 0.03 

19  Mar  - 4 Apr  1975 

199 

03-21  h 

- 

1.1  f 0 . 1 

The  results  obtained  here  are  listed  in  Table  5.  The  values  of  A are 
given  for  each  section  of  Fig  11  and  for  Fig  12,  together  with  an  estimate  of  the 
height  and  local  time  at  which  the  values  apply.  The  'seasonal  bias',  as  given 
by  the  latitude  of  perigee  and  time  of  year  (see  Fig  11),  is  also  indicated.  The 
height  is  taken  as  above  the  average  perigee  height"*  . 

The  first  section  of  the  fitted  curve  in  Fig  II  gives  a value  of  A = 1.1, 
and  during  nearly  the  whole  of  this  period  the  perigee  is  experiencing  'winter' 
conditions.  For  the  second  section,  with  A = 0.9,  the  perigee  is  enjoying  pre- 
dominantly 'summer'  conditions;  for  the  third  and  fourth  sections,  where  A = 1.0, 
perigee  has  average  seasonal  conditions.  (Although  the  third  section  may  seem 
to  have  a bias  towards  'winter'  the  short  'summer'  occurs  during  the  time  of 
greatest  change  in  i .)  These  results  suggest  that  A may  be  higher  in  winter 
than  in  summer.  Such  a seasonal  dependence  of  A has  not  been  detected  before 
in  the  analysis  of  satellite  orbits,  but  measurements  at  specific  sites  by  the 
radar  back-scatter  method  ' 1 ’ *b ' ' ' indicate  that  the  strongest  west-to-east  winds 
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usually  occur  in  the  winter,  is  higher  in  the  winter, 

lend  support  to  this  conclusion. 


The  results  here 
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The  first  four  values  of  in  Table  5 (for  average  local  time)  are  all 

below  the  average  curve  in  Tig  I 1 of  Ket  33.  This  decrease  in  the  value  of  A 
can  also  be  detected  in  other  analyses  of  satellite  orbits*^  *l  between  1971  and 
1976,  when  compared  with  the  average  curve  of  Ref  33,  which  is  based  mainly  on 
results  from  the  1960s.  So  it  seems  possible  that  the  average  rotation  rate  of 
the  upper  atmosphere  has  decreased  during  the  early  1970s:  further  evaluations 
of  atmospheric  rotation  rate  are  needed  before  this  tentative  conclusion  can  be 
confirmed. 

The  value  obtained  for  the  last  15  days  in  orbit,  A = 1.1  + 0.1  at  a height 
of  199  km,  is  at  a time  close  to  equinox,  but  it  is  not  averaged  over  all  local 
times  like  the  other  values.  The  local  time  runs  from  03  h through  midnight  to 
2!  h and  there  is  a mild  bias  towards  the  evening  hours  because  the  decrease  in 
inclination  is  then  more  rapid.  So  the  value  might  be  described  as  'average  with 
a slight  bias  towards  evening',  and  indicates  west-to-east  winds  of  25  t 25  m/s 
at  an  average  latitude  of  58°N.  Results  for  such  high  latitudes  are  unusual  in 
satellite  orbital  analysis,  because  the  effect  of  atmospheric  rotation  becomes 
very  small  as  perigee  approaches  apex  (maximum  latitude)  and  it  is  not  usually 
possible  to  obtain  accurate  values  of  A . The  successful  results  obtained  here 
are  attributable  to  (a)  the  large  change  in  orbital  period  in  the  last  few  days 
before  decay,  and  (b)  the  frequent  and  accurate  orbits  obtained  with  the  aid  of 
the  NORAD  observations. 

7 DKNS1TY  SCALK  HEIGHT 

7. I Introduction 

The  perigee  distance,  a(l  - e)  , gradually  decreases  under  the  influence  ol 

air  drag,  and  the  decrease  is  proportional  to  the  density  scale  height,  11  . So 

values  of  H can  be  found  from  the  decrease  in  a(l  - e)  . If  the  odd  zonal 

16 

harmonic  and  lunisolar  perturbations  are  removed  from  a(l  - e)  , using  TROD  , 
the  remaining  variation  should  show  a steady  decrease  as  a result  of  air  drag 
alone.  Those  ' correc ted ' values  of  perigee  distance,  Q , for  the  85  orbits  of 
Table  2 are  plotted  as  crosses  in  Fig  3 of  Ref  7,  with  a smooth  curve  drawn 
through  the  points.  This  graph  has  been  reproduced  here  as  Fig  13.  The  values 
of  a(l  - e)  for  the  15  daily  orbits  before  decay  (Table  4)  are  plotted  as 
089  circles  in  Fig  14,  and  after  removal  of  per turbat ions  are  plotted  as  crosses  with 
a smooth  curve  drawn  through  the  values.  The  curve  of  Fig  14  can  be  regarded  as 
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superseding  that  of  Fig  13  from  MJD  42492  onwards.  The  variation  of  Q , in 
Figs  13  and  14,  shows  a fairly  smooth  decrease  due  to  the  action  of  air  drag. 

7 _ 2 Me t hod 


The  theoretical  equation  for  the  variation  of  Q is 


2h.m2  i 


7 I - 2e  + 


sin  i cos 


, n / 2 H HJ  e 1 \ ) 

u + >a>  2 2’  2 J j * 


where  H is  the  density  scale  height,  H is  the  value  of  H at  perigee,  11^ 

is  the  value  of  H at  a height  1.5  H above  perigee,  and  e'  is  the  ellipticity 

of  the  atmosphere,  taken  as  equal  to  the  Earth's  ellipticity,  0.00335. 

Equation  (34)  is  valid  for  ae/H  > 3. 

P 

Average  values  of  Q were  calculated  over  a time-interval.  At  , long 
enough  to  ensure  that  an  accurately  measurable  change,  A Q , in  Q , has  occurred. 
The  values  of  AQ/At  serve  as  values  of  Q in  equation  (34),  and  averaged  values 
of  Mj  , M,,  e,  a,  i and  cos  2m  over  the  corresponding  At  were  used  to  cal- 
culate values  of  H ^ from  equation  (34). 

For  1971-106A,  the  value  of  ae/H  exceeds  3 from  launch  until 

P 

MJD  42495,  and  17  values  of  were  calculated  during  this  time:  they  are 

plotted  as  circles  in  Fig  15,  and  the  time  over  which  they  are  averaged  is 
indicated.  The  corresponding  average  values  of  height  Vj(  = y + '-5  1U  are 
given  at  the  top  of  Fig  15.  The  standard  deviations  in  the  values  of  derive 

from  the  estimated  errors  in  the  values  of  AQ  arising  from  errors  in  Q . The 
errors  due  to  neglected  terms  in  134)  are  smaller  than  the  errors  in  AQ  , and 
are  ignored. 


Alter  M.1D  -»2495,  when  the  value  of  ae/H  falls  below  3,  equation  (.34) 

P 42,43 

becomes  inaccurate,  and  the  small-eccentricity  form  of  the  theory  " ’ must  be 
used.  The  theoretical  equation  for  Hj  when  ae/H^  <3  is  H|  = AQ/f(z)  , 
where  z " ae/H.  and 


sph .atm. 


Values  of  ^r^  - r^j^H  for  a spherical  atmosphere  are  given  in  Fig  18  of  Ref  42; 

and  jijifz)  - iMZj)^  . the  oblateness  correction,  is  given  in  Fig  22  of  Ref  42.  The 
values  of  H^  apply  at  a height  ■yH^  above  perigee,  where  > is  given  in  Fig  9 
of  Ref  43.  Seven  values  of  H were  evaluated  between  MJD  42495  and  decay,  the 
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last  four  being  daily  values.  For  these  seven  values,  the  observed  values  of  e 

must  be  corrected  for  the  odd-harmonic  oscillation,  and  this  was  done  using  PROD^  3 . 

These  seven  values  of  H are  also  plotted  as  circles  in  Fig  15  with  their 

standard  deviations.  The  sd  is  the  estimated  error  in  the  value  of  AQ  (between 

1.2%  and  3.8%)  arising  from  errors  in  Q , together  with  the  estimated  error  in 

reading  Figs  18  and  22  of  Ref  42  (between  2%  and  6%).  The  errors  due  to  neglected 

42 

terms  in  equation  (35)  should  be  less  than  1%,  and  are  ignored.  The  correspond- 
ing average  values  of  height  yj(=  Yp  + YHp)  are  given  at  the  top  of  Fig  15. 


The  values  of  H.  plotted  as  crosses  in  Fig  15  are  the  values  of  density 
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scale  height  obtained  from  the  COSPAR  International  Reference  Atmosphere  1972 
for  heights  y^  and  the  appropriate  exospheric  temperatures,  , obtained 

from  Ref  7. 


7 . 3 Discussion  of  results 

The  general  impression  gained  from  studying  Fig  15  is  that  the  ClRA  1972 
values  are  in  fairly  good  agreement  with  the  values  from  the  orbital  analysis. 

This  confirms  the  accepted  opinion  that  ClRA  1972  represents  the  scale  height  H , 
a very  variable  parameter,  fairly  accurately.  The  average  difference  between  the 
observational  and  ClRA  1972  values  in  Fig  15  is  about  1%,  while  the  rms  differ- 
ence is  about  10%. 


Three  of  the  observational  values  in  Fig  15  differ  from  the  ClRA  values  by 
more  than  5 km,  those  at  MJD  41859-42018  (1973  June  26  - December  2),  MJD  42253- 
42329  (1974  July  25  - October  9)  and  MJD  42398-42438  (1974  December  17  - 
1975  January  26). 

Tlie  first  of  these  values  is  low  compared  with  that  given  by  ClRA,  probably 
because  the  density  is  low  over  this  period,  relative  to  other  years  (see  Fig  14 
of  Ref  7).  The  second  value,  which  is  higher  than  ClRA,  is  at  a time  when  density 
is  fairly  average  (see  Fig  14  of  Ref  7),  though  it  is  rising  towards  the  end  of 
the  period;  so  no  reason  for  this  high  value  can  be  given.  The  third  value, 
from  1974  December  17  to  1975  January  26,  which  is  also  higher  than  ClRA,  is 
however  explained  by  the  above-average  values  of  density  over  the  whole  period 
(see  Fig  8 of  Ref  7). 

The  results  also  show  that  it  is  possible  to  obtain  accurate  and  consistent 
daily  values  of  H at  the  end  of  the  satellite's  life  when  good  orbits  are 
available  from  numerous  observations.  These  last  few  values  are  also  in  good 
agreement  with  ClRA  1972. 
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8 CONCLUSIONS 

The  orbit  of  197I-I06A  has  been  determined  at  85  epochs  during  its  40-month 
life  from  6635  observations.  The  accurate  Hewitt  camera  observations  were  used 
in  24  of  these  orbits,  in  which  the  average  sd  in  inclination  corresponds  to 
about  80  m in  distance;  while  for  the  other  61  orbits  the  average  sd  corresponds 
to  about  170  m in  distance.  The  other  orbital  parameters  show  improved  accuracy 
with  the  camera  observations,  but  to  a lesser  degree.  The  residuals  of  the 
observations  have  been  sent  to  the  observers. 

A further  2000  observations  were  received  from  NORAD  covering  the  last 
15  days  of  the  satellite's  life.  With  these  observations,  15  daily  orbits  were 
determined.  These  orbits  at  the  end  of  the  satellite's  life  are  more  accurate 
and  at  closer  intervals  than  has  been  possible  in  the  past;  the  average  sd  in 
inclination  corresponds  to  about  120  ra  in  distance. 

Analysis  of  the  inclination  and  eccentricity  at  14:1,  29:2  and  15:1 
resonances,  has  yielded  lumped  values  of  the  harmonic  coefficients  of  order  14, 

29  and  15.  The  lumped  14th-order  coefficients  from  fitting  i and  e with 
SIMKES  are  given  in  equation  (16)  and  have  been  used  in  a recent  evaluation  of 
the  14th-order  harmonics.  The  analysis  of  the  change  in  inclination  at  29:2 
resonance  yielded  lumped  coefficients,  equation  (22),  which  should  be  useful 
in  the  future  in  determining  the  individual  coefficients  of  order  29.  The 
variation  in  eccentr ic i ty  was  too  small  to  analyse  successfully.  The  S1MRES 
fitting  of  inclination  and  eccentricity  at  15th-order  resonance  was  very  satis- 
factory. The  lumped  coefficients  are  given  in  equation  (31)  and  will  be  useful 
in  a further  determination  of  the  individual  coefficients  of  15th  order. 

The  variation  of  inclination  between  the  resonances  has  been  analysed  to 
obtain  four  values  of  the  average  atmospheric  rotation  rate,  A . The  results 
are  summarized  in  Table  5.  Because  of  the  very  slow  movement  of  perigee,  two  of 
the  values  of  A had  strong  seasonal  bias  and  the  results  suggest  that  A may 
be  higher  in  winter  than  in  summer.  This  seasonal  difference  has  not  been  detec- 
ted before  from  analysis  of  satellite  orbits,  but  is  in  agreement  with  measure- 
ments at  specific  sites  by  the  radar  back-scatter  method.  The  values  of  A are 
lower  than  expected,  and  it  seems  possible  that  the  average  rotation  rate  of  the 
upper  atmosphere  was  lower  during  the  early  1970s  than  in  the  1960s.  The  daily 
orbits  over  the  last  15  days  of  the  satellite's  life  give  west-to-east  winds  of 
25  + 25  m/ s at  an  average  latitude  of  58°N.  Results  for  such  high  latitude  are 
unusual,  but  are  possible  here  due  to  the  large  change  in  orbital  period  in  the 
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last  tew  days  before  decay,  and  the  frequent  and  accurate  orbits  obtained  with 
the  aid  of  NORAD  observations. 

Values  of  density  scale  height,  Hj  have  been  determined  from  analysis  of 
the  variation  of  perigee  height.  A total  of  24  values  was  obtained  during  the 
satellite's  40-month  life:  eight  of  these  were  in  the  last  15  days  of  the  life 
and  the  last  four  were  daily.  The  values  obtained  show  fairly  good  agreement 
with  CIRA  1972;  the  average  difference  between  the  observational  and  CIRA  1972 
values  is  about  1%  and  the  rms  difference  1 0% . Only  three  observational  values 
differ  from  CIRA  1 9 .7’  by  more  than  5 km  and  reasons  for  two  of  these  differences 
have  been  suggested.  The  agreement  shows  that  CIRA  1972  provides  a good  measure 
of  the  density  scale  height  for  the  years  1972-1975. 
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Table  2 (continued) 
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Table  2 (concluded) 
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Fig  5 Variation  of  inclination  at  14th-order  resonance 
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Fig  7 Variation  of  inclination  at  29:2  resonance 
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Fig  9 Variation  of  inclination  at  15th-order  resonance 
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Variation  of  eccentricity  at  15th-order  resonance 


Fig  1 1 Values  of  inclination  cleared  of  perturbations 
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Values  of  a(1  — e)  for  Cosmos  462 


